IN THE CLAIMS : 

Please amend the claims to read as follows: 
1.-20. (Canceled) 

2 1 . (Previously presented) An adaptive direct volume rendering system, 
comprising: 

one or more central processing units for executing programs; 
a user interface for receiving a plurality of volume rendering parameters; and 
an adaptive volume rendering engine module executable by the one or more central 
processing units, the module comprising: 

instructions for fragmenting a sampled 3-D dataset of a scalar field into a 
plurality of sub-volumes of different sizes, each sub-volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent upon 
data values of the scalar field and an illumination model; 

instructions for selectively casting a plurality of rays from a 2-D image plane 
towards the sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 

instructions for selecting a subset of the plurality of sub-volumes for 
interacting with the ray; 

instructions for estimating the ray energy reflected by each sub-volume 
of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel value at 
the selected location on the 2-D image plane; and 

instructions for estimating pixel values at other locations on the 2-D 
image plane using the pixel values at the selected locations 

wherein the instructions for selectively casting a plurality of rays from a 2-D 
image plane include 

subdividing the 2-D image plane into a plurality of sub-planes; 
for each of the plurality of sub-planes, 

casting four rays from the four corners of the sub-plane and estimating 
a pixel value at each corner; 
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calculating a maximum pixel value variation within the sub-plane; and 
recursively subdividing the sub-plane into multiple child sub-planes of 
smaller sizes by casting a ray from the center of the sub-plane until the maximum pixel value 
variation of the sub-plane is below a predefined imaging error threshold. 

22. (Original) The system of claim 2 1 , wherein 

the step of fragmenting a sampled 3-D dataset includes 

fragmenting the 3-D dataset into eight sub- volumes; and 
for each sub-volume, recursively fragmenting it into eight smaller sub- 
volumes until the size of the smallest sub-volumes reaches a predefined size limit. 

23 . (Original) The system of claim 22, wherein the predefined size limit is a sub- 
volume comprising 2x2x2 3-D cells and the eight corners of each cell are associated with 
eight data values of the scalar field. 

24. (Original) The system of claim 23, wherein the data value at any location 
within the cell is tri-linearly interpolated using the eight data values at the eight corners of the 
cell. 

25. (Original) The system of claim 2 1 , wherein the set of parameters include a 
maximum, an average value, and a minimum data value of the scalar field within the sub- 
volume. 

26. (Original) The system of claim 2 1 , further comprises 
instructions for constructing an octree comprising a root node, a plurality of 

intermediate nodes, and a plurality of leaf nodes; 

instructions for associating the root node with the 3-D dataset; 

instructions for associating each of the plurality of leaf nodes with a smallest sub- 
volume from the plurality of sub-volumes; and 

instructions for associating each of the plurality of intermediate nodes with a sub- 
volume from the plurality of sub-volumes that is larger than the smallest sub-volume. 

27. (Canceled) 

28. (Previously presented) The system of claim 2 1 , wherein the maximum pixel 
value variation within the sub-plane is defined as the maximum deviation of pixel values at 
the four corners of the sub-plane from an average pixel value of the sub-plane. 

29. (Previously presented) The system of claim 2 1 , wherein the predefined 
imaging error threshold is modulated by an image rendering speed provided by a user, a 
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distance to an edge of an object embedded in the 3-D dataset, and a difference between a 
pixel value estimated from an adaptive ray casting and a pixel value estimated from a bi- 
linear interpolation. 

30. (Previously presented) An adaptive direct volume rendering system, 
comprising: 

one or more central processing units for executing programs; 
a user interface for receiving a plurality of volume rendering parameters; and 
an adaptive volume rendering engine module executable by the one or more 
central processing units, the module comprising: 

instructions for fragmenting a sampled 3-D dataset of a scalar field into 
a plurality of sub-volumes of different sizes, each sub-volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent 
upon data values of the scalar field and an illumination model; 

instructions for selectively casting a plurality of rays from a 2-D image 
plane towards the sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 
instructions for selecting a subset of the plurality of sub- 
volumes for interacting with the ray; 

instructions for estimating the ray energy reflected by each sub- 
volume of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel 
value at the selected location on the 2-D image plane; and 

instructions for estimating pixel values at other locations on the 2-D 
image plane using the pixel values at the selected locations 

wherein the instructions for selecting a subset of the plurality of sub-volumes for 
interacting with the ray include 

identifying a largest sub-volume along the ray path and its corresponding 
maximum and minimum data values; 

checking if the opacity transfer function varies monotonically between the 
maximum and minimum scalar field values; 
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if the function does not vary monotonically, recursively 

identifying a smaller sub-volume along the ray path and its 
corresponding maximum and minimum data values; and 

checking if the opacity transfer function varies monotonically between 
the maximum and minimum scalar field values of the smaller sub-volume; and 

if the function does vary monotonically, calculating the amount of ray energy 
reflected by the sub-volume during its interaction with the ray. 

3 1 . (Original) The system of claim 30, wherein two lookup tables are constructed 
for the opacity transfer function such that a forward lookup table contains the data value 
difference to a nearest local extreme of the opacity transfer function along the data value 
increasing direction and a backward lookup table contains the data value difference to a 
nearest local extreme of the opacity transfer function along the data value decreasing 
direction. 

32. (Original) The system of claim 3 1 , wherein if the maximum data value of the 
sub-volume is smaller than the summation of the minimum data value of the sub-volume and 
its corresponding data value difference stored in the forward lookup table or the minimum 
data value of the sub-volume is larger than the difference between the maximum data value 
of the sub-volume and its corresponding data value difference stored in the backward lookup 
table, the opacity transfer function varies monotonically between the minimum and maximum 
data values. 

33. (Previously presented) An adaptive direct volume rendering system, 
comprising: 

one or more central processing units for executing programs; 
a user interface for receiving a plurality of volume rendering parameters; and 
an adaptive volume rendering engine module executable by the one or more 
central processing units, the module comprising: 

instructions for fragmenting a sampled 3-D dataset of a scalar field into 
a plurality of sub-volumes of different sizes, each sub-volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent 
upon data values of the scalar field and an illumination model; 
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instructions for selectively casting a plurality of rays from a 2-D image 
plane towards the sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 
instructions for selecting a subset of the plurality of sub- 
volumes for interacting with the ray; 

instructions for estimating the ray energy reflected by each sub- 
volume of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel 
value at the selected location on the 2-D image plane; and 

instructions for estimating pixel values at other locations on the 2-D 
image plane using the pixel values at the selected locations; 

wherein the instructions for estimating the ray energy reflected by each sub-volume of 
the subset include: 

estimating a maximum energy differential of the sub-volume; 

comparing the maximum energy differential against a predefined energy error 

threshold; 

if the maximum energy differential is above the predefined energy error 
threshold, recursively; 

selecting a smaller sub- volume along the ray path; and 
estimating a new maximum energy differential of the smaller sub- 
volume; and 

if the maximum energy differential is below the predefined energy error 
threshold, calculating the amount of ray energy reflected by the sub-volume using the 
illumination model. 

34. (Original) The system of claim 33, wherein the maximum energy differential 
depends on the opacity transfer function and the maximum, average, and minimum data 
values of the sub-volume. 

35. (Original) The system of claim 33, wherein the amount of ray energy reflected 
by the sub-volume depends on the length of ray path within the sub-volume, the opacity 
transfer function within the sub-volume, the average scalar field value of the sub-volume, and 
the local gradient vector of scalar field within the sub-volume. 
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36. (Original) The system of claim 33, wherein if the sub-volume is a smallest 
sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3-D cell within the 
smallest sub-volume, and if the opacity function does not vary monotonically within the cell, 
the 3-D cell is further divided into multiple sub-cells until the dimension of a smallest sub- 
cell reaches the cross-section of the ray. 

37. (Original) The system of claim 33, wherein if the sub-volume is a smallest 
sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3-D cell within the 
smallest sub-volume, and if the opacity function varies monotonically within the cell, the 
maximum energy differential of the 3-D cell is estimated by dividing the maximum energy 
differential of the sub-volume by 2. 

38. (Original) The system of claim 33, wherein if the sub-volume is a smallest 
sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3-D cell within the 
smallest sub-volume, and if the opacity transfer function varies monotonically within the cell 
and an iso-surface exists in the 3-D cell, the maximum energy differential of the 3-D cell is 
calculated using the eight data values at the corners of the 3-D cell and the opacity transfer 
function. 

39. (Original) The system of claim 33, wherein the predefined energy error 
threshold is modulated by an image rendering speed specified by a user, and a zoom factor in 
the case of parallel projection or a perspective angle and a perspective distance between the 
image plane and the 3-D dataset in the case of perspective projection. 

40. (Original) The system of claim 2 1 , wherein 

the instructions for estimating pixel values at other locations on the 2-D image 

plane include 

for each location, 

selecting four pixel values associated with four ray origins surrounding 

the location; and 

bi-linearly interpolating a pixel value at the location using the four 

pixel values. 

4 1 . (Previously presented) A computer program product for use in conjunction 
with a computer system, the computer program product comprising a computer readable 
storage medium and a computer program mechanism embedded therein, the computer 
program mechanism comprising: 
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instructions for fragmenting a sampled 3-D dataset of a scalar field into a 
plurality of sub-volumes of different sizes, each sub-volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent upon 
data values of the scalar field and an illumination model; 

instructions for selectively casting a plurality of rays from a 2-D image plane 
towards the sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 

instructions for selecting a subset of the plurality of sub-volumes for 
interacting with the ray; 

instructions for estimating the ray energy reflected by each sub-volume 
of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel value at 
the selected location on the 2-D image plane; and 

instructions for estimating pixel values at other locations on the 2-D image plane 
using the pixel values at the selected locations 

wherein the instructions for selectively casting a plurality of rays from a 2-D 
image plane include 

subdividing the 2-D image plane into a plurality of sub-planes; 
for each of the plurality of sub-planes, 

casting four rays from the four corners of the sub-plane and estimating 
a pixel value at each corner; 

calculating a maximum pixel value variation within the sub-plane; and 
recursively subdividing the sub-plane into multiple child sub-planes of 
smaller sizes by casting a ray from the center of the sub-plane until the maximum pixel value 
variation of the sub-plane is below a predefined imaging error threshold. 

42. (Original) The computer program product of claim 4 1 , wherein 
the step of fragmenting a sampled 3-D dataset includes 

fragmenting the 3-D dataset into eight sub-volumes; and 
for each sub-volume, recursively fragmenting it into eight smaller sub- 
volumes until the size of the smallest sub-volumes reaches a predefined size limit. 
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43. (Original) The computer program product of claim 42, wherein the predefined 
size limit is a sub-volume comprising 2x2x2 3-D cells and the eight corners of each cell are 
associated with eight data values of the scalar field. 

44. (Original) The computer program product of claim 43, wherein the data value 
at any location within the cell is tri-linearly interpolated using the eight data values at the 
eight corners of the cell. 

45. (Original) The computer program product of claim 41, wherein the set of 
parameters include a maximum, an average value, and a minimum data value of the scalar 
field within the sub-volume. 

46. (Original) The computer program product of claim 41 , further comprises 
instructions for constructing an octree comprising a root node, a plurality of 

intermediate nodes, and a plurality of leaf nodes; 

instructions for associating the root node with the 3-D dataset; 

instructions for associating each of the plurality of leaf nodes with a smallest sub- 
volume from the plurality of sub-volumes; and 

instructions for associating each of the plurality of intermediate nodes with a sub- 
volume from the plurality of sub-volumes that is larger than the smallest sub-volume. 

47. (Canceled) 

48. (Previously presented) The computer program product of claim 41 , wherein 
the maximum pixel value variation within the sub-plane is defined as the maximum deviation 
of pixel values at the four corners of the sub-plane from an average pixel value of the sub- 
plane. 

49. (Previously presented) The computer program product of claim 41 , wherein 
the predefined imaging error threshold is modulated by an image rendering speed provided 
by a user, a distance to an edge of an object embedded in the 3-D dataset, and a difference 
between a pixel value estimated from an adaptive ray casting and a pixel value estimated 
from a bi-linear interpolation. 

50. (Previously presented) A computer program product for use in conjunction 
with a computer system, the computer program product comprising a computer readable 
storage medium and a computer program mechanism embedded therein, the computer 
program mechanism comprising: 
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instructions for fragmenting a sampled 3-D dataset of a scalar field into a 
plurality of sub-volumes of different sizes, each sub-volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent upon 
data values of the scalar field and an illumination model; 

instructions for selectively casting a plurality of rays from a 2-D image plane 
towards the sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 

instructions for selecting a subset of the plurality of sub-volumes for 
interacting with the ray; 

instructions for estimating the ray energy reflected by each sub-volume 
of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel value at 
the selected location on the 2-D image plane; and 

instructions for estimating pixel values at other locations on the 2-D image 
plane using the pixel values at the selected locations 

wherein the instructions for selecting a subset of the plurality of sub-volumes 
for interacting with the ray include 

identifying a largest sub-volume along the ray path and its corresponding 
maximum and minimum data values; 

checking if the opacity transfer function varies monotonically between the 
maximum and minimum scalar field values; 

if the function does not vary monotonically, recursively 

identifying a smaller sub-volume along the ray path and its 
corresponding maximum and minimum data values; and 

checking if the opacity transfer function varies monotonically between 
the maximum and minimum scalar field values of the smaller sub-volume; and 

if the function does vary monotonically, calculating the amount of ray energy 
reflected by the sub-volume during its interaction with the ray. 

5 1 . (Original) The computer program product of claim 50, wherein two lookup 
tables are constructed for the opacity transfer function such that a forward lookup table 
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contains the data value difference to a nearest local extreme of the opacity transfer function 
along the data value increasing direction and a backward lookup table contains the data value 
difference to a nearest local extreme of the opacity transfer function along the data value 
decreasing direction. 

52. (Original) The computer program product of claim 5 1 , wherein if the 
maximum data value of the sub-volume is smaller than the summation of the minimum data 
value of the sub-volume and its corresponding data value difference stored in the forward 
lookup table or the minimum data value of the sub-volume is larger than the difference 
between the maximum data value of the sub-volume and its corresponding data value 
difference stored in the backward lookup table, the opacity transfer function varies 
monotonically between the minimum and maximum data values. 

53. (Previously presented) A computer program product for use in conjunction 
with a computer system, the computer program product comprising a computer readable 
storage medium and a computer program mechanism embedded therein, the computer 
program mechanism comprising: 

instructions for fragmenting a sampled 3-D dataset of a scalar field into a 
plurality of sub-volumes of different sizes, each sub-volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent upon 
data values of the scalar field and an illumination model; 

instructions for selectively casting a plurality of rays from a 2-D image plane 
towards the sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 

instructions for selecting a subset of the plurality of sub-volumes for 
interacting with the ray; 

instructions for estimating the ray energy reflected by each sub-volume 
of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel value at 
the selected location on the 2-D image plane; and 

instructions for estimating pixel values at other locations on the 2-D image 
plane using the pixel values at the selected locations 
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wherein the instructions for estimating the ray energy reflected by each sub- 
volume of the subset include 

estimating a maximum energy differential of the sub-volume; 

comparing the maximum energy differential against a predefined energy error 

threshold; 

if the maximum energy differential is above the predefined energy error 
threshold, recursively 

selecting a smaller sub-volume along the ray path; and 
estimating a new maximum energy differential of the smaller sub- 
volume; and 

if the maximum energy differential is below the predefined energy error 
threshold, calculating the amount of ray energy reflected by the sub-volume using the 
illumination model. 

54. (Original) The computer program product of claim 53, wherein the maximum 
energy differential depends on the opacity transfer function and the maximum, average, and 
minimum data values of the sub-volume. 

55. (Original) The computer program product of claim 53, wherein the amount of 
ray energy reflected by the sub-volume depends on the length of ray path within the sub- 
volume, the opacity transfer function within the sub-volume, the average scalar field value of 
the sub-volume, and the local gradient vector of scalar field within the sub-volume. 

56. (Original) The computer program product of claim 53, wherein if the sub- 
volume is a smallest sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3- 
D cell within the smallest sub-volume, and if the opacity function does not vary 
monotonically within the cell, the 3-D cell is further divided into multiple sub-cells until the 
dimension of a smallest sub-cell reaches the cross-section of the ray. 

57. (Original) The computer program product of claim 53, wherein if the sub- 
volume is a smallest sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3- 
D cell within the smallest sub-volume, and if the opacity function varies monotonically 
within the cell, the maximum energy differential of the 3-D cell is estimated by dividing the 
maximum energy differential of the sub- volume by 2. 

58. (Original) The computer program product of claim 53, wherein if the sub- 
volume is a smallest sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3- 
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D cell within the smallest sub-volume, and if the opacity transfer function varies 
monotonically within the cell and an iso-surface exists in the 3-D cell, the maximum energy 
differential of the 3-D cell is calculated using the eight data values at the corners of the 3-D 
cell and the opacity transfer function. 

59. (Original) The computer program product of claim 53, wherein the predefined 
energy error threshold is modulated by an image rendering speed specified by a user, and a 
zoom factor in the case of parallel projection or a perspective angle and a perspective 
distance between the image plane and the 3-D dataset in the case of perspective projection. 

60. (Original) The computer program product of claim 41 , wherein 

the instructions for estimating pixel values at other locations on the 2-D image plane 

include 

for each location, 

selecting four pixel values associated with four ray origins surrounding the 

location; and 

bi-linearly interpolating a pixel value at the location using the four pixel 

values. 

61.-68. (Canceled) 
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